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This paper summarizes the status of the finite-element approach to gauge theories, in which the Heisenberg 
operator equations of motion are converted to operator difference equations consistent with unitarity. A review of the 
CO 'entire program, from quantum mechanics to quantum field theory, is given in Ref. pi. 
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ABSTRACT 

We apply the finite-element lattice equations of motion for quantum electrodynamics 
to an examination of anomalies in the current operators. By taking explicit lattice 
divergences of the vector and axial- vector currents we compute the vector and axial- 
vector anomalies in two and four dimensions. We examine anomalous commutators 
of the currents to compute divergent and finite Schwinger terms. And, using free 
lattice propagators, we compute the vacuum polarization in two dimensions and 
hence the anomaly in the Schwinger model. 
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We begin by reminding the reader of the form of the free finite-clement lattice Dirac equation: 
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- ^~W'm,n+l — ' l hs,n) + "^"('/ , m j + l,rni,i 

- VWj.m.L.n) - Mfen = 0. (1) 



. Here \x is the electron mass, h is the temporal lattice spacing, A is the spatial lattice spacing, m represents a spatial 

dattice coordinate, n a temporal coordinate, and overbars signify forward averaging: 
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It is a straightforward exercise to show that, apart from a contact term, the free electron Green's function may be 
expressed as 
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Here the lattice momentum is given by 
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p = — , Q = p° = \/p 2 +l^ 2 , (tp)i = t&npiTr/M. (4) 



The quantity O is related to Co by 
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c = -tan— . (5) 
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Figure 1: Plot of the lattice vacuum polarization Il(0)/(e 2 /7r) for M = 2533 as a function of v = //A/2. Shown are 
curves with r = h/A = 0.1, 0.75, 0.99. 



We take M , the number of lattice points in a given spatial direction, to be odd, so that ip is periodic on the spatial 
lattice. 

A similar expression can be derived for the free photon propagator. 

This propagator (||) has been used to perform a calculation of the vacuum polarization in two-dimensional QED, 
with a result consistent with the anomaly in the Schwinger model, e 2 /ir. Some typical results are shown in Fig. 1. 
A similar calculation of the anomaly in four-dimensional electrodynamics is in progress. 

Interactions 

Interactions of an electron with a background electromagnetic field is given in terms of a transfer matrix T: 

ip n+ i = T n ip n , (6) 
which is to be understood as a matrix equation in m. Explicitly, in the gauge ^4° = 0, 
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where 

< m < = (-l)^ + ^[e m „ S co S C m ,, S . 

+ ^sinC^r^JsecC^ <W,m^- (8) 

Here 

!1, to' > m, 
0, to' = m, (9) 
— 1, TO' < TO, 

and (the following are local and unaveraged in m^, n) 

A M 

^ = f 4tt ; c w) =E^' do) 



and 
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Cm^mr = Sgn ( m i ~ TO j) S S n ( m j ~ m jKm'J> (H) 
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with 

sgn (to - m') = e m . m , - 8 m , m >- (12) 

Because T> is anti-Hermitian, it follows that T is unitary, that is, that <f> m ,n = ^Am,n is the canonical field variable 
satisfying the canonical anticommutation relations. 



It is instructive, and very simple, to consider the Schwinger model, that is the case with dimension d = 2 and 
mass fi — 0. We set h — A because the light-cone aligns with the lattice in that case. Then we see that the transfer 
matrix for positive or negative chirality, that is, eigenvalue of 175 = 7°7 1 equal to ±1, is 



T ± = i±£ (13) 



It is an immediate consequence of M) that 



(T+) m , m , = <W+ie 2iCm , (14) 
(T_) m , m / = 5m+i,m'e _2lCm ', (15) 

which simply says that the + (— ) chirality fermions move on the light-cone to the right (left), acquiring a phase 
proportional to the vector potential. 

Using ( |l4| ) and ( |l5| ) the anomaly in the Schwinger model is easily computed. The corresponding calculation in 
four-dimensional electrodynamics will be presented elsewhere. 
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Figure 2: Coefficients of three spectral fit for different lattice sizes, (/i = 0) For the first derivative term, the coefficient 
shown is iSA 2 . 



Anomalous Current Commutators 

It is extremely interesting to compute commutators of the gauge-invariant lattice current, and compare with the 
anomalous commutators in the continuum: 

([J°(0, x), J(0)]> = iSVS{x) + idVV 2 5{x), (16) 

where S is the quadratically divergent Schwinger term [^|, and d = 1/12-7T 2 for the Bjorken-Johnson-Low regularization 
|||. We have carried out a straightfoward evaluation of the current-current commutators on the lattice, and have 
performed a fit with lattice delta functions. The results of such fits for the coefficients of the first three odd derivatives 
of delta functions are shown in Fig. 2 for various lattice sizes. The coefficient of the first derivative is roughly consistent 
with the Schwinger result . Similarly, the result for the coefficient of the third derivative term in the commutator is 
in decent agreement with the BJL result (3); however, in both cases, there seem to be significant discrepancies. 
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